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ABSTRACT 


A  number  of  subprograms  are  read  from  a  magnetic  tape  into  the  storage 
of  an  electronic  computer  during  the  execution  of  a  computing  job.  The 
order  in  which  these  subprograms  are  arranged  on  the  tape  is  being  studied, 

A  number  of  jobs  are  executed  during  a  certain  period  of  time  using  different 
groups  of  subprograms  from  the  tape.  It  is  our  aim  to  re-arrange  the  sub¬ 
programs  on  the  tape  in  order  to  reduce  the  time  it  takes  to  read  the  sub¬ 
programs  into  storage. 
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FOREWORD 


An  expression  for  the  total  time  of  selection  of  the  subprograms  is 
determined*  It  is  not  the  aim  of  this  paper  to  find  the  ideal  distribution 
in  order  to  minimize  the  selection  time*  Instead  a  method  of  trial  and 
error  is  used  to  improve  the  distribution  sufficiently  for  practical  pur- 
poses*  Exchanges  between  two  subprograms  at  a  time  are  made*  If  an  ex¬ 
change  does  not  reduce  the  time  of  selection  a  re-exchange  is  made*  A  pro¬ 
gram  in  FORTRAN  for  the  IBM  7090  is  supplied  to  complete  the  task*  A 
practical  example  is  given. 
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SELECTION  OF  SUBPROGRAMS 


A  subprogram  on  magnetic  tape  consists  of  a  number  of  records  with 
record  gaps  between  them.  Each  record  contains  a  certain  number  of  "words". 
A  word  has  a  fixed  length  on  the  tape.  Hence,  the  length  of  a  subprogram 
may  be  measured  In  terms  of  the  number  of  words  It  contains  within  each 
record,  allowing  for  record  gaps.  Let  us  suppose  we  have  n  subprograms 
Si»88,...,s^  of  lengths  ti,t2,...,t^  words  on  the  tape.  The  "head",  which 

reads  the  words  from  the  tape.  Is  located  at  the  beginning  of  the  tape,  and 
this  Is  also  the  beginning  of  Sj^ .  To  select  the  subprogram  s^  (m  >  1  )  the 

head  moves  over  (tj^  +  tg  +  ...,  +  t^  ^  )  words.  After  the  first  subprogram 

s^  has  been  read  Into  storage^  the  head  Is  located  at  the  beginning  of  • 

To  select  the  next  subprogram  s^j  the  head  moves  over 

a)  (t  .  +  t  ,  +...,+  t  )  words  If  r  >  m  +  1 

ro+i  mrs  r-l 

b)  (t  +  t  +  t  )  words  if  r  <  m. 

'  '  m  m-i  *  r' 

It  takes  no  time  to  select  s  (r  =  m  +  1).  If  the  same  subprogram  is 

used  again  (r  =  m)  ve  can  ignore  the  corresponding  selection  time,  since 
It  cannot  be  reduced  by  changing  the  distribution.  In  this  fashion^ the 
head  moves  forwards  and  backwards  over  the  tape  to  select  all  the  subpro¬ 
grams  for  a  job. 

We  assume  that  the  head  moves  forwards  and  backwards  over  the  tape 
with  the  same  constant  speed.  Hence,  the  total  selection  time  may  be 
measured  in  terms  of  the  number  of  words  passed  by  the  head  without  read¬ 
ing. 
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II 


TOTAL  TIME  OF  SELECTION 


There  are  n  subprograms  Sj  (i  on  the  tape  of  lengths 

tj  (i  =  l,...,n).  During  a  certain  period  of  tlrae^  a  number  of  jobs  have 
been  executed  and  the  numbers  f^]((l,k  »  and  fj  (1  = 

have  been  recorded,  fj  Is  the  number  of  times  the  subprogram  Sj  was 
selected  first  during  the  execution  of  a  job.  fjjc  Is  the  number  of  times 
the  selection  of  s^  was  followed  by  the  selection  of  ,  To  find  the 
total  time  of  selection,  T,  we  form  the  following  expressions: 

fi(ti  +  t^  +  ...,  +  1  >  1 

+  i  +  ^+3  tk-2  ^ 


^kl(^  •••>  ^ 


and  find  the  sum  of  all  possible  terms: 

n  1-1  *^*2  n 

T  =  }  ft  )  ;  ftk 

1  s  3  m=i  1=1  k  si^a 


k-i 


-I  1  I  I  I  I 

«  =  j.  ms7  i  ^  n  ir  ^  o  ms|  irso  isi  m^i 


m 


k  =  2  1=1 


We  change  this  express loi*  for  faster  computation.  The  third  term  may 
be  changed  to 


n-i 


^  n  n-i  k  -1  k 

m=t 


t  . 

m 


k  s  8  1=2 

The  second  term  may  be  changed  to 


k=3  1=2 


m=i 


n-2 


n-3  n 


1  =  1 


k  •! 


L  ""LX  ^**'1 


1=  1  k=i+3 


m 
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The  first  term  may  be  written  as 


tx  + 


n-x 

I 


t+1 


1 


mssi 


m 


Hence , 

n-i 

T  =  fg  •  tj  +  2^  (fljl  + 

k  =  2  »*! 


m=si 


n-1 


n  n-l 


1  +  1  H 


i  =  2 


m=i 


1  =  2 


n-l  k 


k  =  3  1  S  3 


n-l  n-l  k 

fl-2  ^k  +  1 

1.1  j «  3  Its  £  tn=l-i 


The  lasttuo  terms  add  up  to 
n-l  n-x 


^  ^  ,1-1  ■*"^i-3,k  +  x)  ^ 


t  . 
m 


1=3  k=  1 


m=  i-x 


This  form  of  '£  has  been  used  In  the  subroutine  FINITE  to  compute  T. 
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REDUCTION  OF  THE  TIME  OF  SELECTION 


To  reduce  the  time  of  selection  T  we  change  the  order  in  which  the 
subprograms  are  on  the  tape.  We  compute  at  first  for  the  original  dis¬ 
tribution.  Then  we  exchange  two  subprograms  Sj  and  Sjj  and  compute  T^  .  If 
Ti  ^  Tq  ,  we  exchange  Sj  and  Sj,  again.  Otherwise  we  make  the  next  exchange. 

We  let  i  =  1,2 . n-1  and  k  =  i  +  1,  i  +  2,...,n  until  a  permanent  change 

has  been  made.  Then  we  start  again  with  i  =  1  and  k  =  2.  This  method  will 

lead  to  a  distribution  which  cannot  be  improved  by  the  exchange  of  any  two 

subprograms.  For  n  >  50  it  may  take  hours  of  computer  time  before  the  final 
distribution  has  been  found,  because  the  changes  become  more  and  more 
insignificant  and  harder  to  find  as  time  goes  on.  It  is  advisable  to  run 
the  program  for  10  minutes  at  a  time  and  break  it  off  as  soon  as  the 
improvement  becomes  insignificant.  In  an  example  with  n  =  65  the  time  T 

was  reduced  by  26%  during  the  first  10  minutes  of  computer  time,  by  7% 

during  the  next  30  minutes  run,  by  4%  during  the  next  30  minutes  run  and 
by  1%  during  the  next  20  minutes  run.  In  an  example  with  n  =  71  the  time 
T  was  reduced  by  5.5%  in  7  minutes  and  by  1,6%  in  another  30  minutes.  In 
this  case  the  program  worked  on  an  improved  distribution  of  subprograms. 

To  enable  a  restart  of  the  program  without  a  change  of  the  input  data 

fji,  ,  f j  ,  and  tj  ,  the  program  transforms  the  matrices  fn  ,  fj  ,  and  tj  to  any 

given  order  of  distribution  s  ,s,  ,s  . s  ,s  ,s  .  We  simply  type  the 

^  d  DCxyz 

integers  a,  b,  y,  z  on  certain  data  cards,  and  the  program 

establishes  this  order  and  continues  making  improvements. 
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IV 


A  PROGRAM  IN  FORTRAN 


A  program  in  FORTRAN  for  the  IBM  7090  is  provided.  The  arrays  F(N,N), 
T(N)  and  S(N)  contain  the  matrices  f^j  >  and  fj  ,  where  N  =  n.  The  array 

^nJ(N)  contains  the  integers  a,  b,  c,...,x,  y,  z  mentioned  in  3.  The  array 
NU(N)  is  initially  made  to  agree  with  M0(N)  and  then  subjected  to  changes, 
indicating  the  current  order  of  distribution^ 

®NU(1)’  ®NU(2)’-"»  ®NU(N)- 


The  total  time  of  selection  is  denoted  by  TT,  the  change  in  time  due 
to  an  exchange  of  subprograms  by  DT.  The  following  subroutines  are  usedj 


READ; 


EXCH: 


HEAD; 

ORDER; 

INIT; 

FINDT; 

PLAN; 


This  subroutine  reads  from  the  data  cards  the  arrays  T,F,MIJ 
and  S.  F  must  be  typed  columnwise  on  the  cards,  e.q.  fxi >^2 
^ni .•••  .  where  f^j  =  f^g  =  ...,  f^^  =  0. 

MU  must  contain  the  numbers  I,  2,  N  at  the  start  and 

later  the  improved  order  of  distribution. 


This  subroutine  exchanges  the  subprograms  S  and  S. 


Lb* 


It 


makes  the  necessary  changes  in  the  matrices  F,T,NU  and  S. 

This  subroutine  prints  the  total  time  TT  and  the  order  of  dis¬ 
tribution  NU  as  a  heading  for  each  page. 

This  subroutine  prints  the  total  time  TT  and  the  order  of  dis¬ 
tribution  NU  whenever  10  changes  have  been  made. 

This  subroutine  transforms  to  the  order  of  distribution  NU  =  MU. 
This  subroutine  evaluates  TT. 

This  subroutine  exchanges  two  subprograms  at  a  time  and  prints 
for  each  permanent  change  the  numbers  Li ,  L2  and  DT»  indicating 
that  the  Li  th  and  L2nd  subprograms  have  been  exchanged  and  that 
TT  has  been  increased  by  DT« 


The  main  program  calls  INIT, FLAN ^ORDER  and  EXIT. 

N  is  taken  to  be  71  in  the  following  listing  of  the  program. 
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ORDER 


MAIN  PROGRAM 

1  CALL  INIT 

2  CALL  PLAN 

3  CALL  ORDER 

4  CALL  EXIT 

END(1 ,0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0) 

SUBROUTINE  INIT 

DDENSION  F(71  ,71  ) ,T(71  )  ,NU(71  ) ,MU(71  ) ,S(71  ) ,LU(71  ) 
COJMON  F,T,TT,DT,NU,LN,N,L1 ,L2,MU,S 

1  N=71 

2  CALL  READ 

17  DO  19  1=1  ,N 

18  NU(I)=I 

19  LU(I)=I 

20  N1=N-1 

21  DO  27  1=1  ,N1 

22  LI =I 

23  L2=MU(I) 

24  L2=:LU(L2) 

25  IF(L2-L1 >24,27,26 

26  CALL  EXCH 

27  LU(I)=L2 

28  CALL  FINDT 
31  CALL  HEAD 

RETURN 

END(1 ,0,0, 0,0, 0,1 ,0,0, 0,0, 0,0, 0,0) 
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SUBROUTINE  EXCH 

DIMENSION  F(71 ,71) ,T(71  ) ,NU(71 ) ,MU(71 ) ,S(71 ) 
COJMON  F,T,TT,DT,NU,LN,N,L1 ,L2,MU,S 

5  D06  1=1  ,N 
A=F(L1 ,1) 

F(L1 ,I)-F(L2,I) 

6  F(L2,I)=A 

7  DOlO  1=1  ,N 

8  A=F(I,L1) 

9  F(I,L1)=F(I,L2) 

10  FU,L2)=A 

1 1  A=T(L1 ) 

12  T(L1)=T(L2) 

13  T(L2)=;A 

14  A=NU(L1) 

15  NU(L1 )=NU(L2) 

16  NU(L2)^ 

17  A=S(L1) 

18  S(L1)=S(L2) 

19  S(L2)^ 

RETURN 

END(1 ,0,0, 0,0, 0,1 ,0,0, 0,0, 0,0, 0,0) 

SUBROUTINE  READ 

DIMENSION  F(71  ,71  ) ,T(71  ) ,NU(71 ) ,MU(71 ) ,S(71  ) 
CXllMON  F,T,TT,DT,NU,LN,N,L1  ,L2,MU,S 

1  READ  INPUT  TAPE  5,2,T 

2  P0RMAT(16F5.0) 

3  READ  INPUT  TAPE  5,4,F 

4  FORMAT(24F3.0/24F3.0/23F3.0) 

5  READ  INPUT  TAPE  5, 6, MU 

6  F0RMAT(40I2) 

7  READ  INPUT  TAPE  5,4,S 
RETURN 

END(1 ,0,0, 0,0, 0,1 ,0,0, 0,0, 0,0, 0,0) 
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SUBROUTINE  FINDT 

DIMENSION  F(71 ,71 ),T(71 ),NU(71 ),MU(71 ),S(71 ) 

COMMON  F,T,  TT,DT,NU,LN,N,L1 ,L2,MU,S 

1  N1=N-1 

2  TT=0.0 

3  A=T(1 ) 

TT=TT+A*S(2) 

4  DO  6  1=2, N1 

5  A=A4T(I) 

6  TT=TT+A*(  F(  1 , 1  )-«  ( 1+1  )  ) 

A=A+r(N) 

TT=TT+A*F(N,1 ) 

7  DO  9  1=2, N1 

8  A=;A-T(I-1 ) 

9  TT=TT+A*F(N,I) 

10  DO  14  1=3, N1 

11  A=T(I-1) 

12  DO  14  J=1  ,N1 
1  3  A=5A+T(  J) 

14  TT=TT+A*(F(I-2,J+1 )+F(J,I-l )) 

15  DO  16  1=2, N1 

1 6  TT=TT+r( I )*F(I- 1  , 1+l  ) 

17  RETURN 

END(1 ,0,0, 0,0, 0,1 ,0,0, 0,0, 0,0, 0,0) 

SUBROUTINE  HEAD 

DIMENSION  F(71 ,71 ),T(71 ),NU(71 ),MU(71 ),S(71 ) 

COMMON  F,T,TT,DT,NU,LN,N,L1 ,L2,MU,S 

1  WRITE  OUTPUT  TAPE  6,2,TT,NU 

2  F0RMAT(1H1 ,9X,5HTIME=, El 6.8,1 2HTHE  ORDER  IS/(1X,I2,39I3)) 

3  LN=0 

4  RETURN 

END(1 ,0,0, 0,0, 0,1 ,0,0, 0,0, 0,0, 0,0) 
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SUBROUTINE  PLAN 

DIMENSION  F(71 ,71 ),T(71 ),NU(71 ),MU(71 ),S(71 ) 

COJtlON  F,T,TT,DT,NU,LN,N,L1 ,L2,MU,S 
A=TT 

1  N2=N-1 

2  DO  16  1=1 ,N2 
L1=I 

3  L3=L1+1 

4  DO  16  J=L3,N 
L2=J 

5  CALL  EXCH 
CALL  FINDT 
DT=TT-A 

6  IF(DT)7,16,16 

7  WRITE  OUTPUT  TAPE  6,9,L1,L2,DT 

9  FORMAT (1  OX, lOHCHANGE  L1  =  ,I4,8H  AND  L2=,I4,4H  DT=,E16.8) 
A=TT 

1 0  LN=LN+1 

11  IF(50-LN)12,12,13 

12  CALL  HEAD 

13  IF(XM0DF(LN,10))2,14,2 

14  CALL  ORDER 

15  GO  TO  2 

16  CALL  EXCH 

17  RETURN 

END(1 ,0,0, 0,0, 0,1 ,0,0, 0,0, 0,0, 0,0) 

SUBROUTINE  ORDER 

DIMENSION  F(71 ,71 ) ,T(71 ) ,NU(71 ) ,MU(71 ),S(71 ) 

COJWON  F,T,TT,DT,NU,LN,N,L1  ,L2,MU,S 

1  WRITE  OUTPUT  TAPE  6,2,TT,NU 

2  PORMAT(8X,13HNEW  0RDER.TT=,E16.8/(1X,I2,39I3)) 

LN=LN+1 

RETURN 

END(1 ,0,0, 0,0, 0,1 ,0,0, 0,0, 0,0, 0,0) 
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